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Abstract 

We have examined the influence of rotation on the potential energy and the transport 
coefficients of the collective motion (friction and mass coefficients) . For axially symmetric 
deformation of nucleus Th we found that at excitations corresponding to temperatures 
T > 1 MeV the shell correction to the liquid drop energy practically does not depend 
on the angular rotation. The friction and mass coefficients obtained within the linear 
response theory for the same nucleus at temperatures larger than 2 MeV are rather stable 
with respect to rotation provided that the contributions from spurious states arising due 
to the violation of rotation symmetry are removed. At smaller excitations both friction 
and mass parameters corresponding to the elongation mode are growing functions of 
rotational frequency u> ro t. 

Keywords: collective motion, rotating nuclei, linear response theory, transport coefficients 
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1 Introduction 



The recent success of Flerov Laboratory, JINR, Dubna in the synthesis of the superheavy 
compound systems with Z = 114, 116 has provoked a considerable theoretical interest to the 
fusion-fission reactions at low excitation energies. Commonly such reactions are described by 
solving the Langevin equation [|l|]-@ for the collective variables which specify the shape of the 
nuclear compound system formed in the result of fusion of heavy ions with nuclei. Usually 
such systems are formed with rather high angular momentum. The effect of rotation on the 
fusion or fission probability is included at most in the calculation of the macroscopic part of 
the deformation energy. The possible dependence on rotation of the shell correction as well as 
friction and inertia is completely ignored. However one might expect the strong dependence 
on rotation of the transport coefficients since the rotation changes considerably the single- 
particle spectrum. To the best of out knowledge this effect was analyzed only in Q where 
rather strong dependence of friction coefficient on angular velocity was found. 

In the present paper we have continued the investigations along the line presented in ||. 
Besides the friction coefficient we have paid also attention to the rotational dependence of the 
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mass parameter and the shell correction. The computations are carried out with two-center 
shell model [||, || which allows for rather flexible parametrization of the shape around the 
touching point and which was used earlier in dynamical computations 0. As the compound 
nucleus we chose the system 208 Pb + 16 O =^> 224 Th for which the experimental information on 
the temperature dependence of the damping parameter is available || and which was studied 
in [[|, [1(J without account of the rotation. 

Due to technical reasons we had to limit ourselves to the excitations above T = lMeV. The 
point is that the rotation violates not only the axial symmetry but the time reversal symmetry 
too and the BCS approximation to pairing interaction breaks down. To account for the pairing 
accurately one has to solve a kind of Hartree-Fock-Bogolyubov equation what requires rather 
time consuming computation. For this reason we considered here the excitations corresponding 
to temperature larger than T = lMeV when the pairing can be neglected. 

The paper is organized as follows: In Section |2] we quote the main relations of the linear 
response theory adopted for the rotating nuclei. The quasi-static properties (moment of 
inertia, liquid drop and shell component of potential energy) are examined in Section [3[ 
The influence of rotation on the response functions and transport coefficients is investigated 
in Secs.[|,|5]. The special attention here is paid to the elimination of the contribution from 
spurious states caused by rotation. The main conclusions and open questions are formulated 
in Summary. The expressions for the matrix elements of J x on two-center oscillator basis wave 
functions are given in the Appendix. 



2 General formalism 

By describing the rotating nuclei one usually transforms the Hamiltonian from the laboratory 
co-ordinate system to the body fixed (or intrinsic) co-ordinate system. In the result, instead 
of the Hamiltonian H{Q^) one has to consider the Routhian operator 

R(Qfi,^rot) = H{Q^) — LdrotJx (1) 

with uj ro t being the rotational frequency and J x - the projection of angular momentum on the 
rotation ). The variables in (|l[) are the deformation parameters which specify 

the shape of the (deformed) mean field. 

The energy of rotating nucleus E = (u rot \H\u rot ) and angular momentum I = (u> rot \J x \uj rat ) 
are growing function of the rotational frequency u rot . At small values of u rot one can use the 
perturbation theory to obtain 

E = Eq + -Jcran^rot ( 2 ) 

where J cr an is the cranking model moment of inertia. In the approximation of independent 
particles it is given by 

J C ran = tfY. T ^ ± \{k \ J X \j)? (3) 
kj e i tk 

The rik in @ are the (temperature dependent) occupation numbers and summation is carried 
out over the single-particle states k and j. If u rot is not small the single-particle spectrum 
and eigen-f unctions are to be found numerically by solving the eigen-values problem, 

R{Q^uJrot)\k) = e k {Q^u rot )\k) (4) 
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Like in the case without rotation the transport coefficient of collective motion, the tensors 
of stiffness C^ u , friction 7^ and mass M^ v can be derived within the linear response theory 



□I O from the so called collective response function Xcoiifw) approximating X co u(u;) by the 



response function of damped oscillator 

[xcoiiHW — ► M-muj 2 - + cy'ku (5) 

Here k is the coupling tensor , see ([□]) below. The collective response function Xcoai^) * s 
related to the Fourier transform X^( u ) °f the intrinsic (causal) response function 

X^(t) = 0(f)£tr (p qs (Q„T)[F^t),Fl(0)]) (6) 

by 

XcquM = *(« + x{")) x xH (7) 

The k in (|7|) is the inverse of coupling tensor ([TT]) and operators i 7 ^ (t) in @ are the interaction 
representation for the derivatives of the Routhian R{Q^,uj rot ) with respect to deformation (or 
rotational frequency u> rot ), 

» [) ' dQ, ' ^ du rot 1 j 

The average in ([]) is calculated with the quasi-static density operator for which the canonical 
distribution is assumed, p qs (Q^,uj rot ,T) oc exp(—R(Q^,u> ro t)/T). 

The response function (H) can be used to calculate the deviation of the average value of 
Ffj, from its quasi-static value (calculated at some deformation point Q°), (see |TT| ) 

POO 

HF») t = -£ / xAt - s)(Q»(s) - Ql)ds (9) 

JJ J —00 

The Fourier transform of @ reads 

«(4 = -EU*(w) (10) 

V 

with 5Q v {uj) being the Fourier transform of (Qj,(s) — 
The coupling tensor fc (c.f.[0]) is 

- (V 1 )^ = - Kvv = c^(o) + x^(o) (11) 

where stiffness C M!/ (0) of the free energy JF(Q,T) and static response x^i®) are defined by 
the static properties of the system, 

C (0) ~ ^ mT) (12) 
and X M z,(0) is the Fourier transform of the intrinsic response function (Q) taken at u = 0. 
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3 Quasi-static properties 



The collective potential energy E(Qp, I) is one of the most essential ingredients appearing 
in the theory of large scale collective motion. The derivatives of E{Q^ L) I) with respect to 
deformation define the collective conservative forces. The second derivatives of E(Q /Ji , I) 



(stiffness) is used to find the inverse of the coupling tensor k^ u (|TT1) which appears in the 
collective response function (|7|). 

Like in the case without rotation we will use for calculation of the potential energy the 
Strutinsky shell correction method [13, The idea of applying the Strutinsky renormal- 



ization to the rotational problem was advanced by Pashkevich et al |~5], jl6| , see also [[TjJ. 



Following [16j one can express the intrinsic energy E(Q fl , I) as 

E(Qp, I) = E LDM {Q„ I) + 5R{Q^ I) (13) 

where E LDM (Q^, I) is the liquid drop energy of rotating nucleus and 5i?(Q M ,J) is the shell 
correction, 



6R = Y,W -V- -/ n[ — )eg{e)de (14) 




In the case of finite temperature instead of the shell correction to the intrinsic energy one has 
to consider the shell correction to the free energy SR ==>- 8T = 5R — T5S, where SS is the shell 
correction to the entropy, see fl9|]. The energies eu in ( |14|) are to be found by diagonalization 
of the shape dependent Routhian ([|) and g(e) is the average density of single-particle states, 
see 0,0]. 

As argued in fllEj] , Eldm(Qh, -0 can De represented rather accurately by 

E LDM (Q^ I) = E LDM (QJ + I 2 /2J rig (QJ (15) 

Here, J r i g {Qp) is the rigid body moment of inertia for the rotation around x-axes and 
Eldm{Qv) is the liquid drop energy of non-rotating nucleus. 

In the computations presented below we will use the two-center shell model |5|, || and 
consider only axially symmetric shapes. Such shapes are specified mainly by two parameters: 
the distance z between the centers of left and right oscillator potentials and the parameter 
S which fix the spheroidal deformation of the "fragments". Below we will consider the case 
when deformations of left and right fragments are the same, Si = 62 = 5. Furthermore we 
will consider here only a one-dimensional path in the deformation space and define 5 = S(zq) 
looking for the minimum of the total energy at fixed zq, see @. 

The comparison of the rigid body and cranking model moment of inertia is shown in 
Fig.l. It is seen that both methods give rather close results (the pairing was neglected). The 
dependence of J cra n on the rotation is also not very strong. So, the main source of rotational 
dependence of E LDM (Q^, I) is the I 2 - term in (|i~5|). 

The left-hand-side part of Fig. 2 shows the rotational dependence of the liquid drop part 
of deformation energy (we suppose that spherical shape corresponds to = 0) 

El e L(Q„ I) = E LDM (Q„ I) - Eldm(Qh = 0, /) (16) 

As it is seen the rotational dependence of the deformation energy is rather strong. The fission 
barrier becomes lower due to rotation and disappears completely at / « 60h for the nucleus 
224 Th shown in the figure. 
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Figure 1: The deformation dependence of the rigid-body J r % g . (curve with dots) and the 
cranking model J cra n moments of inertia. The deformation parameter zq is here the distance 
between the centers of left and right oscillator potentials of two-center shell model and m is 
the nucleon mass. The Jcran is computed for the temperature T = 1 MeV. 



The effect of rotation on the fission barriers is known for decades and taken into account 
nowadays in all computations of the deformation energy. The rotational dependence of the 
shell correction is less clear. The diagonalization of Routhian ([!]) is much more time consuming 
due to the break of axial symmetry as compared with the diagonalization of the non-rotating 
shell model Hamiltonian. Hence it is assumed usually that this dependence is weak and the 
shell correction is computed at uj rot = only. To clarify this point we have computed the shell 
correction for several values of J as a function of deformation along the liquid drop fission 
valley of 224 Th. Indeed, see right-hand-side of Fig. 2 , the fluctuation of 6 J 7 is less then 1 MeV 
for variation of / from zero to I = 60h. Very likely such weak dependence of 8T on / can be 
neglected. 

The weak dependence of the shell correction on rotation is not so surprising. It was 



pointed out by Strutinsky [18] that the shell effects are not sensitive to rotation as far as 
the perturbation hu rot is small compared with the spacing hfl between the gross shells, 
Tiuj rot <C hQ , with hQ ^8 — 10 MeV. For the maximal value of spin / = 60h shown in Fig. 2 
the fiuJ ro t = hi I J ~ 0.4 MeV at the saddle and hu rot « 0.6 MeV at the minimum. Both 
values of hu ro t are much smaller than HQq. Since the moment of inertia J increases with 
growing deformation the huj ro t gets smaller (for fixed /). This explains, at least partly, why 
at large deformation the shell correction is less sensitive to rotation as compared with small 
deformation, see Fig.2. 



4 Response functions 

The intrinsic response function x M! ,(cg>) is one of the most simple and important quantity of 
the linear response theory. The friction and mass coefficients in the so called zero frequency 
approximation are expressed in terms of derivatives of the intrinsic response functions, see 
ITTJ . The intrinsic response function is also an important ingredient of the collective response 
function (|^). So we will look first at the effect of rotation on the intrinsic response function. 
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Figure 2: The liquid drop deformation energy (left) and the shell correction 8T = 5R — TSS to 
the free energy (right) for temperature T = 1 MeV as function of the deformation parameter 
zq. The solid, dashed, dotted-dashed and dotted lines correspond to the values of angular 
momentum equal to 0, 20, 40 and 60ft. 



4.1 Intrinsic response function 

The Fourier transform of the intrinsic response function given by @ can be expressed as the 
sum over single-particle states 

xM = J2x 1 k(^)F^ J (17) 

jk 



with 



— - n(Q) (Qk{£l)Gj(£l + uJ + ie) + gj(Q)g k (Q - u - 



is 



Here n(Q) is the Fermi function determining the occupation of the (rotation dependent) single- 
particle levels. The Qy. appearing in QT8| ) is the one-body Green function 

g k (uj±ie) = - — 1 (19) 

Too — e k — h'(u, T) ± il(u, T)/2 

It is parameterized in terms of the real and imaginary part of the self-energy E(cu,T) = 
H'(u,T) — iT(u,T)/2. The T(u,T) is assumed to have the form 

r(W ' T) ~ T Q 1 + [{hu - W + ^] /c 2 (20) 

and H'(u,T) is coupled to T(u,T) by the Kramers-Kronig relation. The Qk(oj) represents the 
distribution of single-particle strength over more complicated states. It is related to Q k by 

Q k (u) = i{Gk{uJ + ie) - Q k {u} - is)) (21) 



For the simplified case when the collisional damping could be neglected the intrinsic re- 
sponse function attains the form 



^ h(uj - uJki) + lO 

with TiuJki = £fc — Q and iO being infinitely small imaginary constant. 

Fig. 3 shows few examples of the response function (|T7|)- (|2T|) for the elongation mode. 
The deformation parameter in this case is the distance q = q(zo, S) between the left and right 
centers of mass (divided by the diameter of the sphere with the same volume) and F g -operator 
is , see M 

p = dH_ dzo dH_d5__(dq_ 86_ dqy 1 (dB_ d6_dH_\ 

q dz dq 85 dq \ dz Q dz d5 J \ dz dz 85 J 

where the derivative dS/dz is to be taken along the fission path 5 = S(z ). 

Comparing the response function corresponding to different values of uj rot one notices a 
peak in the low frequency region. This peak is absent in the case u) rot = 0. With growing 
uj rot the peak gets "stronger" and moves away from uj = 0. Its position is approximately 
proportional to to rot . This circumstance hints that this additional peak may be caused by 
rotation (let us call it here "rotational" peak). Recalling that the friction coefficient 7 and 
mass parameter M (at least in the zero frequency limit) are defined by the derivatives of 
response function with respect to uj at uj = it is clear that the value of both 7 and M can 
be very sensitive to the "rotational" peak. The numerical results show that the contribution 
from "rotational" peak to the real part of the response function can lead to a negative value 
of mass parameter. Thus, the "rotational" peak could be of spurious origin and one has to 
treat this problem very accurately in order to get the reliable results for friction and inertia. 



4.2 The conservation of angular momentum 

We will consider in this section the case of a single deformation parameter Q for simplicity. 
The generalization to multi-dimensional case is straightforward. 

It is clear that the Routhian R(Q, uJ ro t) = H(Q) —uJ ro tJx violates the rotational symmetry. 
The operator J of angular momentum does not commute with Routhian, thus J is not a good 
quantum number. By varying uj rot one can fix the average value of J x , i.e. one chooses uj rot 
in such a way that 

(uj rot \J x \uj rot ) = I (24) 

where by \uJ rot ) we denote the Irast state of Routhian (|1|) calculated for given uj TOt . 
If in addition to rotation we switch on also the vibrations 

H(Q) - uj rot J x H(Q°) - uj rot J x + F5Q(t) (25) 

then \uJrot) becomes time-dependent, \uJ rot ) =^ \^rot) t and, in principle, average value of J x is 
not conserved any more 

(uJrot\Jx\Urot) t = I(t) = 1 + 8{J x ) t 1 (26) 

The variation S(J x ) t is given within the linear response theory by (H), namely 

rod 

5(J x ) t = - / x JmF (t - s)(Q(s) - Q°)ds (27) 
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Figure 3: The frequency dependence of the real (left) and imaginary (right) parts of F q F q 
intrinsic response function for several values of the rotational frequency uo rot . The dotted, 
dashed and solid curves correspond to temperatures T = 1, 2 and 3 MeV. The heavy 
solid curve marks the modified response function Xff(uS), see (p5|), for T = 1 MeV. The 
computations are done for the "ground state" shape of nucleus 22A Th very close to the sphere, 
z = 0.1, parameters of spheroidal deformations 5\ = 62 = 0.05. 
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The 5{J x ) t given by fl2"T|) is not zero. The possible way to make {u rot \ J x \oj TOt ) t time independent 
is to allow u rot to dependent on time 

L0 rot UJ rat (t) = UJ° rot + 8u rat (t) 

H{Q)-u r J x H(Q°)-co° rot J x + F5Q(t)-JJco rot (t) (28) 

The time- dependent correction Soo ro t(t) should be found from the requirement that S(J x ) t (or 
its Fourier transform ^(Jx)^) is equal to zero, 

5(J x ) t = 6(J X )„ = (29) 

This problem can be easily solved by means of linear response theory. Considering the time- 
dependent part F5Q{t) — J x 5uj rot {t) as a small perturbation one can find (mind @) 

8(J*)t = - / XJ x F(t-s)(Q(s)-Q°)ds 

Xj x J x (t - s)(uw(s) - uj° rot )ds, (30) 



or its Fourier transform 

s (Jx) u = -Xj x f{u)5Q{u) - XJ^M^roti^) (31) 
The analogous expression can be also written for the variation S{F) uj 

S(F) W = -xff(u)5Q(uj) - XFjM5u ro t{u) (32) 
Recalling (p9|) one can find 8uj rot {uj) from (|3lD as 

8u) rot {u) = -Xj x f{u)8Q(u)/xj x jM (33) 

In principle, 5u rot (t) could be found by Fourier transform of (|33|). But this is not necessary 
for practical purpose. One can insert 5uj rot {uj) in ([32|) to define the modified response function 

= -Xff(u})5Q(uj) (34) 

with 

Xff(^) = Xff(^) r~\ — ( 35 ) 

XJxjAu) 



The modified response function (|35|) is shown by heavy solid line in Fig. 3. It is seen that 
"rotational peak" has disappeared. Consequently, the transport coefficients computed with 
modified response function fl35|) will not contain the contributions from the "rotational peak" 
and would differ considerably from these derived with xff^)- 

The above method was successfully used in ]19[ to remove the contributions to the response 
function caused by the violation of the particle number conservation by pairing. It was 
demonstrated there that fixing of the average value of the particle number with the time- 
dependent density matrix leads to the same secular equation for the vibrational mode as 
obtained within RPA. This method is rather general and can be used to fix any physical 
quantity. For example, using transformation analogous to (^) one can remove the center of 
mass motion in case of isoscalar dipole vibrations. 



9 



It is of certain interest to compare the secular equation which results from the above 



approach with one obtained earlier, for example within the so called "cranked RPA" [[20] , pi 



It is argued in |20], |2lJ that for the description of non-rotational excitations in the rotating 
nuclei one should substitute the Routhian ([]]) by some supplementary rotational invariant 
Hamiltonian H 

R = H — oj rot J x =>■ H = H — h(J 2 ) (36) 

The many-body Hamiltonian H was approximated in |20], ^TJ by the mean-field part H$ phis 
quadrupole-quadrupole interaction 



m=2 



H = H -- £ (-l) m Q 2m Q 2 - 



(37) 



and for h(J ) the expansion was used 



h(J ) = (h)+ u rot (J x - (J x )) + fi x (J x - (J x )f + \i{Jl + J 2 Z ) + ... 



(3? 



with 



u rot = 2fi(J x ), fi x = ^d 2 (H)/d(J x ) 2 



(39) 



The terms omitted in (|3^) are of the higher order in J x — (J x ), J y or J z . Solving of the 
Hamiltonian H within RPA leads to the secular equation 



with 



to 



uj 2 T + (uj) = 



S X X S X Q 



(40) 



S 



x2 



Sq x «Soo — k/2 Si 



02 



2x 



S- 



20 



S 22 - k/2 



(41) 



where 



S, 



\- f (wfc - ni)u)ig jn) » ( n k ~ ni)vki Jy,) Jv) 

/ j \ , ,2 _ , ,2 "kl Hkl , 2 _ , ,2 "kl Hkl 



(42) 

k>i I w " - w w ~ u= ki ) 

The spurious zero energy state is separated by (|40| ) and all other excitations are given by 
F+(u) = 0. 

To get the equation for the collective frequencies in our approach let us recall that the 
equation (|29"D together with the self-consistency condition 

S(P) U = k5Q(u) (43) 

leads to the system of equations (mind (|32|) , (|3T]) ) 



(k + xff(u))5Q(u) + XFj x (u)5u rot (u) 








(44) 



10 



The eigenfrequencies for the system (H3) are found from the equation 



Det(w) = (k + Xff(u))xj x j x (uj) - XfjMxjM^) 







(45) 



If one would neglect the effects of collisional damping then one could use expressions ( fZ2"| ) for 
the response functions. In this case Det(u;) coincides exactly with the nontrivial part J z+ (uo) 
of the secular equation (|I0D obtained in |20[ (one should also put S% v = = since only 
axially symmetric shapes are considered in present work). So Eq.(|45"D does not contain the 
spurious contributions caused by the violation of rotational symmetry. 
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Figure 4: The friction and inertia (|j) for the elongation mode as functions of the average value 
of angular momentum / (upper x-axes) or rotational frequency u rot (lower x-axes). Dotted, 
dashed and solid curves correspond to temperatures T = 1, 2 and 3 MeV . The computations 
are done for several deformations of nucleus 224 Th which correspond to the minimum of 
potential energy, saddle and two touching spheres. 



5 Transport coefficients 

In Fig. 4 we show the friction ^ qq and mass M qq coefficients defined by the fit @ of the 
collective response function for three particular deformations which are of a special interest: 
at the ground state deformation of 22A Th, at the saddle and for at touching point. The 
last configuration is of the interest for description of initial stage of fusion reactions. The 
deformation parameter q = Ri 2 /2R is here the distance R12 (divided by the diameter of 
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nucleus) between the centers of mass of left and right parts of nucleus. The advantage of 
such a choice is explained in |], |lOl. As it is seen from Fig. 4 the dependence of 7 gg on uj rot 
is much weaker as that found in [ffl. Evidently, this is because we have removed the spurious 



"rotational" peak from the response function. Without modification fl35|) we would obtain the 
friction coefficient which look very much like that of 0]. 

One can see from Fig. 4 that for the ground state deformation and temperatures T = lMeV 
the friction and mass coefficients depend somewhat on the rotation. Both 7 W and M qq increase 
with u rot in the interval < u rot < 0.08o;o. For higher T and u) rot both 7 9(? and M qq are rather 
stable with respect to variation of T and uj rot . For more deformed shapes the friction and 
mass coefficients are not very sensitive to the rotation for all temperatures. 
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Figure 5: The reduced friction coefficient f3 qq = ^qq/Mqq (left) and the damping factor r\ qq 



"fqq/2y\C qq \M qq (right) versus temperature. The dotted, dashed and solid curves correspond 
to the values of angular momentum equal to 0, 40 and 60 h. 



Finally, Fig. 5 shows the reduced friction coefficient (3 qq = ■y qq /M qq and the damping factor 
Vqq = lqq/^\J\Cqq\M qq at the saddle of 22A Th as the function of temperature. The damping 
factor reveals whether collective motion is underdamped (rj < 1) or overdamped (rj > 1). As 
Fig. 5 shows, the collective motion changes from underdamped to overdamped at T « lMeV. 
Both f3 qq and r\ qq shown in Fig.5 increase with the temperature. This behaviour is in a 
qualitative agreement with the one found in ||. The increase of r\ qq with the temperature 
is impossible to explain neither with the wall friction nor with the hydrodynamical viscosity. 
The increase of rj qq with temperature obtained here is not as rapid as that found in ||. The 
account of rotation does not diminish this discrepancy since both f3 qq and rj qq shown in Fig.5 
do not depend much on the rotation. The dependence of r] qq on / seen from Fig.5 is mainly due 
to some dependence of liquid drop stiffness on rotation. We should also note that in this work 
we did not normalize the high temperature limit of the mass parameter to the irrotational flow 
value. That is why the numerical results for the mass parameter and, consequently, reduced 
friction (3 qq and damping factor rj qq differ somewhat from these obtained in [H] . 
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6 Summary and outlook 



We have examined the influence of rotation on the transport coefficients of the collective 
motion. Rather unexpectedly we have found out that friction 7 and mass M parameters for 
rotating nuclei are rather sensitive to such fine effects as the violation of rotational symmetry 
by Coriolis term —u rot J x . For the ground state deformation the spurious contributions to 
collective friction and mass are (at least) as large as those of physical importance. This 
circumstance was not clear (to the best of our knowledge) up to now. 

In order to remove the spurious contributions we had to modify the model of "stationary 
rotation" and to introduce the time- dependent rotational frequency. In the result we obtained 
the friction and the mass which demonstrate rather reasonable dependence on the rotational 
frequency uj rot . For excitations above T = 2MeV when the microscopic shell effects disappear 
both 7 and M are rather insensitive to the rotation, i.e. behave like macroscopic quantities. 
For lMeV < T < 2MeV we found some increase of 7 and M with growing uo rot . Such effect 
might be caused by the change of shell structure due to the re-arrangement of single-particle 
states by rotation. 

Even stronger dependence of 7 and M on uo rot should be expected for T < lMeV when 
both shell end pairing effects are present. As it was shown in flT9| , |22fl the pairing effects 
change considerably the collective transport at low excitations. The destruction of pairing 
by the rotation can have considerable effect on the transport coefficients and is worth to be 
examined. Such details could be very important for the accurate description, for example, of 
the final stage of the fusion reaction and formation of superheavy elements which takes place 
at low excitation energy. 

The extension of the method developed in the present work to the simultaneous treatment 
of both pairing and rotation will be the subject of future work. 
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A The matrix elements of J x operator 

The single-particle basis wave functions are defined in the two-center shell model [^ [| as 

\n z n p m s z ) = p nz ( z ) Vm(<p) Xi/zfa) (46) 

where 

Vm &) = -L e im * , (47) 
y An 

with m being an arbitrary integer number and 

i i m-f- 1 

x!?(p) = (-iH 



0(n \ I |m|+l , „ 



^ (n p + |m|) ! 

Here k p = m uj p /h, n p is a not-negative integer and Lj^'(£) is a Laguerre polynomial. The 
z— components of wave function 

(z ,( N~\ U(-n zl - |, -V2k~i(z - si)) for z < 
Vn z [z) \ AT-x U{ - nz2 - 1, ^2kT 2 (z-z 2 )) for z > 
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are expressed in terms of parabolic cylinder functions U(a,x), 

rT( \ _ iFi(a/2 + l/4,l/2 >a; 2 /2) 2/4 

U W> X ) 2 «/2+i/4 T(a/2 + 3/4) 



y^x 1 ^/2 + 3/4 I 3/2 2 xVg _ x . 2/4 

2 a/ 2 +i/4 T(a/2 + l/4) 6 1 j 

The constants iV nzl and iV nz2 are defined by the normalization and the continuity of (f nz {z) 
and its first derivative at z — 0, see [|J. 

The operator of single-particle angular momentum j is 

j = l + s (51) 

where 1 is the orbital momentum 1 = — i [r V] and s is the spin s = ~ cr, with a being the 



Pauli matrices. From (51) it follows immediately that 

(n z n p m s z \j x \n' z n' p m' s' z ) = (n z n p m \l x \n' z n' p m' )5 S2jS / 2 + {s z \s x \s' z )8 n ^8 np ^ p 8 m ^ (52) 
For the spin part of fl52] ) it is easy to find 

(s z \s x \s' z ) = ~(<f,;, a ,+i + <J«;,,.-i) (53) 

The x-component of the orbital momentum l x is given in the cylindrical co-ordinate system 
{r, 9, <p} by 

r .[( d d \ . zcosy? <9 1 

Note that f* = -l x . 

If the nucleus is not left-right symmetric then the rotation axes does not go through z = 
but through the center of mass z cm . Consequently eq.(^4j) should be modified to 

lx = M I (z- Zcm)^ - J smcp H cosl P~q^ I ( 55 ) 

Since the single particle wave functions are separable in {p, z, <p} the matrix elements 
(n z n p m \l x \n' z n' p m' ) are then the product of one-dimensional matrix elements 

(n z n p m \l x \n' z n' p m' ) = 

d 

(n p m\p\n'm') (m\ — isin^|m }{n z \ — \n z ) + (n z \z — z cm \n' z ) x 



dz 

d Id 

— \n p m)(m\i sin (p\m) + (n p m\ — \n' p m')(m\i cos (p—\ 



(56) 



The matrix elements (n z \z — z cm \n' z ) and (n z \d/dz\n' z ) are the same as computed in the two 
center shell model code || ^j. For (m\ — isunp\m') and (m\i cos ipd/d(p\m') it is easy to find 

(m\ -ismip\m') = ^(<5 m ', m+ i -<W,m-i), 
(m\i cos if— \m) = — — (8 m >, m+1 + 8 m >, m -i) (57) 
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What is left to calculate are "p" -matrix elements. These can be calculated using the recur- 
rence relation between Laguerre polynomials and their derivatives. After somewhat lengthy 
derivation one can find 



[n z n p m \l x \n' z n'm' ) 



Jn p + m 6 m f >m 



d 

k p (n z \z - z cm \n' z ) + —f= ( n *I^K) 



5n' p ,n p -i + 



yjn'p + m! 5 m ', m+ i + ^Jr7 p 6. 



'm'm—l 



k p (n^Z z cm\ n z)) 



1 / I 8 I '\ 



5m> 



N' p ,N p +l 



(5? 



Here we have introduced the quantum number N p = 2n p + \m\. The expression fl58|) is 
valid for m and m' being it both non-negative. The matrix elements for non-positive m and 
m! can be related to ( |58|) using (47) and symmetry properties of 4 

— m \l x \n' z n' p — m! ) — (n z n p m \l x \n' z n' p m! ), (59) 

The operator j x couples the states with Aj z — ±1. In this way the states with positive 
and negative j z are coupled to each other. One can reduce the dimension of matrix to be 
diagonalized by factor two using so called Goodman transformation pj. It was suggested in 



23 1 to introduce the basis states of the type 



K) = ^=(| k)+ | k}), 



K) = ^(| k)- | k)) 



(60) 



where | k) =\ n z n p ms z ) for such m and s z that m + s z — 1/2 = j 2 — 1/2 is even. The single 
particle states | n z n p ms z ) with j z — 1/2 - odd up to a sign factor coincide with | fc). It was 
shown in p3| that the matrix elements (K\j x \K') are zero and the matrix of j x on the states 
| K) is of quasi-diagonal form. The nonzero matrix elements are 



(K\j x \K') = -(K\j x \K') = (k\j x \W) 



(61) 



So the quantities of interest are matrix elements (k\j x \k') = (k\j x T\k'}. It turns out possible 
to express the matrix elements {K\j x \K') in terms of (k\j x \k'), 



(K\j x \K') = (n z n p ms z \j x \n' z n'm's' x ) 



(62) 



For (n z n p ms z \j x \n' z n' p m's' z } the expressions (|52]),(|53| ) ,(|5^ ) are to be used. 

The particular case of m + s z = m! + s' z = 1/2 should be considered separately. In this 
case one has to combine (|58f) with (|59D to obtain 



(n. 



n p , 1/2 - s z , s z \j x T\n' z , n' 1/2 - s' z , s' z ) = -5 n ^ n > z 5 np ^ 5, 



s'l- 



-s' z -s z ON' p ,N p -l\jn p + 1/2 - .s. 



,8n> 



1 d 

'k p (n z \z - z cm \n' z ) + -= (n z \ — \n z ) 

>k, 



+-5 s ^ s J N > ptNp+1 Jn' p + 1/2 - s' z 



k p (n z \z z cm \n ) ■= (n z 

h 
« i, 



dz 

d_ 

dz 



(63) 
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